Nonlinear Decoherence in Quantum State Preparation 

of a Trapped Ion 



O 

o 
o 

(N 



00 



> : 
^t- . 

(N ■ 

cn : 
o . 
o ■ 

o ■ 

^ : 
Oh: 

-I— > ' 
G ■ 
cd 

3 : 



Le-Man Kuang, Hao-Sheng Zeng, and Zhao- Yang Tong 
Department of Physics, Hunan Normal University, Changsha 410081, China 

We present a nonlinear decoherence model which models decoherence effect caused by various 
decohereing sources in a quantum system through a nonlinear coupling between the system and its 
environment, and apply it to investigating decoherence in nonclassical motional states of a single 
trapped ion. We obtain an exactly analytic solution of the model and find very good agreement 
with experimental results for the population decay rate of a single trapped ion observed in the NIST 
experiments by Meekhof and coworkers (D. M. Meekhof, et al, Phys. Rev. Lett. 76, 1796 (1996)). 
PACS numbers: 32.80.Pj, 42.50.Lc, 03.65.Bz, 05.45. +b 



In recent years, much progress has been made in 
preparation, manipulation, and measurement of quan- 
tum states of the center-of-mass vibrational motion of a 
single trapped ion experimentally [1-8] and theoretically 
[9-16], which are not only of fundamental physical in- 
terest but also of practical use for sensitive detection of 
weak signals [17] and quantum computation in an ion 
trap [3,9]. In particular, the NIST group [4] has ex- 
perimentally created and observed nonclassical motional 
states of a single trapped ion. In the NIST experiments 
[4], an anti Jaynes-Cummings model (JCM) interaction 
between the internal and motional states of a trapped ion 
is realized through stimulated Raman transitions, which 
couple internal states of the trapped ion to its motional 
states, when the Lamb-Dicke limit is satisfied and the 
driving laser fields are tuned to the first blue sideband. 
Detection of motional states is carried out by observing 
the evolution characteristics of quantum dynamics of in- 
ternal levels of the trapped ion under influence of the anti 
JCM-typed interaction. The NIST experiments revealed 
the fact that the population of the low atomic state (PjJ 
evolves according to the following phcnomenological ex- 
pression 

n(*) = \i i + e^^w^tiv 7 " 4 } (i) 

n 

where p n is the initial probability distribution of motional 
states of the trapped ion in the Fock representation, g 
is a coupling constant between the atomic internal and 
motional states, 7„ is a decay rate. The experimentally 
observed decay rate is of the following form 



7n = 7o(w + l Y 



(2) 



where the observed value of v is v = 0.7. 

A question that naturally arises is: how to explain the 
above experimentally observed decay rate? It is gener- 
ally accepted that the appearance of the decay factor j n 
in the evolution of internal states is a consequence of 
decoherence. It is of practical significance to well under- 
stand decoherence for preparation of nonclassical states 



and quantum computation in ion traps. There are var- 
ious sources of decoherence [1], such as ion vibrational 
decoherence, ion internal-state decoherence, decoherence 
caused by non-ideal external fields, and so on. Recently, 
Schneider and Milburn [18] have investigated decoher- 
ence due to laser intensity and phase fluctuations and 
obtained the power v in Eq.(2) being v = 0.5 instead 
of the experimentally observed value 0.7. More recently, 
Murao and Knight [19], using master equation method, 
have studied decoherence due to the imperfect dipole 
transitions and fluctuation of vibrational potential in the 
NIST experiments. In spite of these efforts, the prob- 
lem of decoherence in quantum state preparation of a 
trapped ion has been not satisfactorily solved, and its 
character and microscopic origin still call for further at- 
tention. In particular, it should be pointed out that the 
experimentally observed decay rate indicated in (2) is a 
collective effect caused by various decohering sources, not 
by a specific decohering source. Nevertheless, authors in 
rcfs[18,19] investigated the decay rate caused only by a 
specific source of decoherence, not by various sources of 
decoherence. So how to model the experimentally ob- 
served decay rate caused by various decohering sources is 
an interesting subject in quantum state preparation and 
manipulation of a trapped ion. In this paper, we present 
a nonlinear decoherence model to model decoherence ef- 
fects caused by various decohering sources in a quantum 
system. We shall show that our theoretical model can 
well describe the experimentally observed decay rate in 
the NIST experiments [4] . 

We consider a single trapped ion with mass m and 
laser cooled to the Lamb-Dicke limit. Following sym- 
bols Ref.[19], we denote three related internal states and 
motional states of the ion by \i) (i = 0, |,|) and \n) 
(n = 0,1,2,...), respectively. The free Hamiltonian of 
the trapped ion is given by Hq = huj x a + a — hwoi\ J,) (J, 
I — hu>Q2\ t)(T |) where ^01(^02) is the transition fre- 
quency between states | J.)(| |)) and |0), a + (a) is the 
creation (annihilation) operator of the motional states 
with the corresponding frequency u> x . Two driving laser 
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beams with detuning A, wave vector ki(k2) and fre- 
quency u>i(u>2) are used to cause dipole transitions be- 
tween the level | |)(| t>) and |0). 

With the dipole and rotating wave approximations, un- 
der large detuning condition the intermediate level |0) 
can be adiabatically eliminated when the Lamb-Dicke 
limit is met and the driving laser beam is tuned to the 
first blue sideband. Then, in the interaction picture of 
Hq, the effective Hamiltonian of the system has the anti 
JCM-typed form 

Hs = %(a +cr + + clcf-) (3) 

where g is a coupling constant, which depends on the 
coupling strength between internal and motional states 
of the trapped ion and the Lamb-Dicke parameter de- 
fined by 77 = Skx n , where 5k is the wave- vector difference 
of the two Raman beams along x, and xq = yjh/2mu x . 
For simplicity, we set Ti = 1 throughout this paper. 

The Hamiltonian (3) is diagonal in the dressed-state 
representation with the following basis 

\ip(n,i)) = -J=(| l,n) - (-1)*| T,n+l»,i = 1,2 (4) 

1^(0,3)) = I T,0) (5) 

And we have Hs\<p(n, i)) — E ni \ip(n,i)) with eigenvalues 
E ni = {-l) i+1 gVn+l for i = 1,2, and E 03 = 0. 

Before going to our model, let us briefly recall a few 
basic facts about the interaction between a quantum sys- 
tem and its environment. The interaction between the 
system and its environment may create two types of ef- 
fects [20-34]: decoherence and dissipation, which can be 
mathematically described by decaying of the off-diagonal 
and diagonal elements of the reduced density operator of 
the system, respectively. These two effects have been 
paid much attention in various areas, for instance, quan- 
tum measurement [20,25-28], condensed matter physics 
[21-23], quantum computation [29-31], and so on. The 
decoherence effect causes the states of the system con- 
tinuously decohere to approach classical states [20,27]. 
The dissipation effect dissipates energy of the system to 
environment [21-23]. The two effects can be understood 
in terms of Hamiltonian formalism [32-34] . If we assume 
that the total Hamiltonian of the system plus environ- 
ment to be Ht = Hs + Hr + Hi, where Hs and Hr 
are Hamiltonians of the system and environment, respec- 
tively, and Hi is the interaction Hamiltonian between 
them, when the Hamiltonian of the system commutes 
with that of the interaction between the system and envi- 
ronment, i.e., [Hs, Hi] = 0, which means that there is no 
energy transfer between the system and the environment, 
energy of the system is conservative, so that what inter- 
action between the system and environment describes is 
decoherence effect. When [Hs,Hi] ^ 0, there is energy 



transfer between the system and environment, so that 
what interaction between the system and environment 
describes is the dissipation effect. It should be pointed 
out that the decoherence and dissipation happen at dif- 
ferent time scales [29,30]. The dissipation effect occurs at 
the relaxation time T re i , while the decoherence time scale 
Td is much shorter than r re i with the time evolution of a 
quantum system. Hence, we here restrict our attention 
on decoherence effect. 

We now present our model. We use a reservoir consist- 
ing of an infinite set of harmonic oscillators to model the 
environment of the single trapped ion in the NIST exper- 
iments, and assume that in the interaction picture of Hq 
the total Hamiltonian is of the following phcnomcnolog- 
ical form 

H T = H s + ]T w k b\b k + F({6 S }) £ c k {b\ + b k ) 
k k 

+^ 2 ({Os})E4- ( fi ) 

k u k 

Here the first term is the Hamiltonian of the system in 
the interaction picture given by Eq.(3); the second term is 
the Hamiltonian of the reservoir; the third one represents 
the interaction between the system and the reservoir with 
a coupling constant c k , where {Os} is a set of linear op- 
erators of the system or their linear combinations in the 
same picture as that of Hs, F({Os}) is an operator func- 
tion of {Os}- In order to enable what the interaction 
between the system and the reservoir describes in Eq.(6) 
is decoherence not dissipation, we require that the lin- 
ear operator Os commutes with the Hamiltonian of the 
system, i.e., [Os,Hs] — 0. It is well known that the 
decohering process can indeed be considered as a quan- 
tum measurement process. The conventional definition 
of a quantum measurement involves any form of inter- 
action between a quantum object and a classical system. 
Therefore, the interaction function F({Os}) in the model 
(6) can involve any form of interaction between the sys- 
tem and environment. This enables it to model the col- 
lective decohering behavior caused by various decoher- 
ing sources. The concrete form of the function F({Os}) 
may be regarded as an experimentally determined quan- 
tity. The last term in Eq.(6) is a renormalization term, 
which is discussed in Ref.[21]. When F({Os}) is a lin- 
ear and nonlinear function of the linear operator Os, We 
call decoherence described by the interaction between the 
system and the reservoir linear and nonlinear decoher- 
ence, respectively, in the similar sense of the linear and 
nonlinear dissipation implied in Rcf.[21]. In this sense, 
decoherence investigated in Ref.[19] is a kind of linear de- 
coherence. In what follow we shall show that nonlinear 
decoherence can better describe the decay rate in the the 
NIST experiments. 

The Hamiltonian (6) can be exactly solved by making 
use of the following unitary transformation 
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After applying the unitary transformation (7) to the 
total Hamiltonian (6), we get a decoupled Hamiltonian 
H' T — Hs + Ysk^kfik, where n k = b\b k . The density 
operator associated with the decoplucd Hamiltonian is 
given by 

p' T (t) = e- iA ^p' T (0)e i6 '^ (8) 

where p' T (0) = ?7/3t(0)C/ _1 , with pt(0) being the initial 
total density operator. Through a converse transforma- 
tion of (7), it is straightforward to obtain the total den- 
sity operator associated with the original Hamiltonian 
(6) with the following expression 

p T (t) = e-^U^e-^^UpTiO) 

xtZ-V^^W^*. (9) 

We assume that the system and reservoir arc ini- 
tially in thermal equilibrium and uncorrelated, so that 
p T (0) = Ps (0) ® Ar(0), where p s (0) and p R (0) are the 
initial density operator of the system and the reservoir, 
respectively. pr(0) can be expressed as pR = Yl k Pk(0) 
where p k (0) = (1 - e ~ 0Uk )e~^ Ukfik is the density oper- 
ator of the fe-th harmonic oscillator in thermal equilib- 
rium, where (3 = l/fc^T, ks and T being the Boltz- 
mann constant and temperature, respectively. After tak- 
ing the trace over the reservoir, from Eq.(9) we can get 
the reduced density operator of the system, denoted by 
p(t) = trnf>T(t), its matrix elements in the dressed state 
representation are explicitly written as 

P(m'A')(m,i)(t) = P(m' ,i')(m,i){Q)Rm'i'mi{t) 

xe-^'i'miW , (10) 
Here the phase is defined by 

0m ; i' mi(^) — \^m'i' ^mi\: (H) 

and Rm'i'mi 

(t) is a reservoir-dependent quantity given by 

R m >i>mi{t) = ~[[Tr R {D(~a mtk )e- lt ^D(-a mk ) 

k 

xD(-a mHlk )e- itUkAk D(-a mHlk )p k (0)}, (12) 

where a mik = f({O mi })c k /uj k with O ml being an eigen- 
value of the linear operator Os in a dressed state, i.e., 
Os\tp(m,i)) = O mi \ip(m,i)), and D(a) = exp(afe^ - 
a*b k ) is a displacement operator. 

Making use of properties of the displacement operator: 

D{a)D{l3) = D(a + /3) exp[Um(a/3*)}, (13) 
exp(xh k )D(a) exp(~xn k ) — exp(ae x b^ — a*e" x b k ) 1 (14) 

and the following formula [35] 



Tr R [D(a)p k (0)] = cxp[-i|«| 2 coth(^)], (15) 

we find that the reservoir-dependent quantity R m 'i'mi(t) 
can be written as the following factorized form 

R m 'i'mi(t) = e-^'.'J'le-r.vJ^ (i 6 ) 
with the following phase shift and damping factor 

&<t>m>Vmi{t) = [F 2 ({O mV }) ~ F 2 {{O mi })} Q 1 (t) , (17) 

J- m'i'mi (t) = [F{{O mli ,}) - F({O mi })} 2 Q 2 (t). (18) 
Here the two reservoir-dependent functions are given by 

Qx(t)= / dojJ(oj)—Y-sm(wt), (19) 
Jo ^ 

Q 2 (t) = 2 r^JH^M sin 2 (^) cothA, (20) 
Jo u z z 

where we have taken the continuum limit of the reservoir 
modes: J2k ~~ * Io° d^-H^)? where J(uj) is the spectral 
density of the reservoir, c{lo) is the corresponding con- 
tinuum expression for c k . 

We assume that the system is initially in a state 
P(0) = I l)U \®T, n Pn\n){n\. Then, from Eqs.(10)-(16) 
we find that at time t the population of the lower atomic 
state is given by 

Pl(t) = ^{l + ^PnCOS^^W+^^Wle^" 1 " 2 ^, 

n 

(21) 

which indicates that the interaction between the system 
and reservoir induces a phase shift S(f) n i n2 (t) and a damp- 
ing factor T n i n2 (t) in the time evolution of the atomic 
population. 

Taking into account the experimental expression (1), 
we choose the following linear operator and interaction 
function: 

O s = a + a+ +acr_, (22) 
F{{Os}) = 6 2 s d+ \ (23) 

where d is an adjustable parameter to describe the non- 
linearity in the interaction, which reflects the deviation 
degree of the nonlincarity of F({Os}) with respect to 
the linear operator Os- The value of the parameter d 
is determined by the experimental results. With these 
choices, it is easy to find that 

F 2 (O nl )-F 2 (O„ 2 )-0, (24) 
F(O nl )-F(O n2 ) = 2(V^TT) 2d+1 . (25) 
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Then the phase shift in Eq.(21) naturally vanishes, and 
the damping factor becomes 

r„i„ 2 =4(n+l) 2d+1 Q 2 (t). (26) 

So that we can find from Eq.(21) that 

Pi(t) = i{i + ^p„cos(2 3 tV^TT) e - 4( " +1) " Q2(t) } 

n 

(27) 

where v = 2d + 1 and Q2(t) is given by Eq.(20). From 
Eq. (27) we see that the argument of the cosine function 
on the RHS of Eq.(27) does have the same form as that 
in the experimental expression (1). Comparing the the- 
oretical expression (27) with the experimental result (1), 
we find that when the nonlinear deviation d = —0.15, 
the n-dependence of the damping factor in Eq.(27) is 
completely in agreement with that seen in the experi- 
mental expression (1). The final step is to determine 
the time dependence of the damping factor in Eq.(27). 
From Eqs.(19), (20), and (27), we see that all necessary 
information about the effects of the environment is con- 
tained in the spectral density of the reservoir. Eq.(27) 
indicates that the time dependence of the damping fac- 
tor is completely determined by the spectral density of 
the reservoir. The experimental expression (1) requires 
that the time dependence of the damping factor must be 
linear, so that if we choose the spectral density such that 

Q2(t) = j-fot (28) 

where 70 is a characteristic parameter, then we can get an 
expression of Pi(t), which has exactly the same form as 
the experimental result (1). It is possible to find a spec- 
tral density of the reservoir to satisfy the condition (28). 
For instance, for the case of zero temperature, if we take 
the spectral density J(co) — 7o/(27rc 2 (w)), substituting it 
to Eq.(20) we can realize Eq.(28). 

In conclusion, we have present a nonlinear decohcr- 
ence model, and obtained its exactly analytic solution. 
It has been shown that our model can give precisely the 
same expression of the population decay rate of the single 
trapped ion as that observed in the NIST experiments [4] . 
The nonlinear dccohercncc model can describe the NIST 
experiments so well. This indicates that the reservior 
and the nonlinear coupling between the system and the 
reservoir, which we design, properly model the real en- 
vironment of the single trapped ion in the experiments. 
It is worthwhile to emphasize that the nonlincarity in 
the coupling describes a collective contribution of vari- 
ous decohering sources to the decay rate. Hence, what 
the nonlinear decoherence describes is a collective deco- 
herence effect caused by various decohering sources not a 
specific decoherence source. We have noted that authors 
in ref.[19] obtained the decay rate in Eq.(l), but deco- 
herence which they considered is a specific decoherence 



caused by the imperfect dipole transitions and fluctua- 
tion of vibrational potential, so their results can not cover 
the contribution of other decohering sources to the decay 
rate in the NIST experiment. It can be expected that the 
nonlinear decoherence model proposed in the present pa- 
per can describe decoherence behaviors of a wide variety 
of quantum systems. 
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